Generation and stability of discrete gap solitons 
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We analyze stability and generation of discrete gap solitons in weakly coupled 
optical waveguides. We demonstrate how both stable and unstable solitons can be 
observed experimentally in the engineered binary waveguide arrays, and also reveal a 
connection between the gap-soliton instabilities and limitations on the mutual beam 
focusing in periodic photonic structures. 
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Suppression of the diffraction-induced beam spreading 
and generation of spatial optical solitons is usually as- 
sociated with the transverse confinement of a beam in a 
self-induced waveguideA. As the medium refractive index 
increases through the nonlinear response, the beam be- 
comes self-trapped in the effective high-index region due 
to the total internal reflection. In periodically modulated 
photonic structures such as fiber Bragg gratings and ar- 
rays of weakly coupled optical waveguides, nonlinearity- 
induced localization of light is also possible in the band- 
gaps which appear due to the Bragg reflection. The later 
mechanism is responsible for the formation of temporal 
Bragg solitons^ and spatial gap solitons 3 . 

Temporal Bragg solitons are usually analyzed in the 
framework of the coupled-mode theory, where the total 
field is presented as a superposition of nonlinearly cou- 
pled counter-propagating waves which experience Bragg 
reflection. This case corresponds to a narrow gap in the 
transmission spectrum, and it is well studied in the the- 
ory of fiber Bragg gratings 2 . 

However, the coupled- mode theory accounts for an iso- 
lated gap only, and such a simplified description of the 
properties of periodic structures is not valid for analyzing 
deeper modulations of the refractive index such as, for ex- 
ample, those in waveguide arrays and photonic crystals. 
Deeper gratings demonstrate more complicated dynam- 
ics, including instabilities of nonlinear guided modes due 
to multi-gap coupling effects and a resonant interaction 
between different bands 4 . 

In this Letter, we study the stability properties of dis- 
crete gap solitons in engineered binary arrays of optical 
waveguides 5 -, which have many common properties with 
spatial gap solitons observed experimentally. Our anal- 
ysis, based on effective discrete equations, captures the 
main features of discrete solitons, and these results are 
confirmed by numerical simulations. In particular, we 
demonstrate how both stable and unstable gap solitons 
can be observed experimentally, and also reveal a funda- 
mental link between soliton instabilities and limitations 
on the mutual beam focusing. 

We consider the head-on excitation geometry, when 
optical beams are incident at grazing angles on a one- 
dimensional periodic structure of weakly coupled optical 
waveguides. Then, the beam self-action and interaction 



resulting in soliton formation can be described by the 
scalar nonlinear Schrodinger equation, 
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where x and z are the transverse and longitudinal coordi- 
nates, respectively, £(x,z) is the normalized electric field 
envelopes, v(x) = v(x + h) is the normalized refractive 
index profile, h is the spatial period, and medium has a 
self-focusing nonlinear response. 

Recently, it has been suggested that the properties of 
spatial gap solitons can be engineered in binary waveg- 
uide arrays, where the first gap is controlled by the width 
difference between alternating "thick" and "thin" waveg- 
uides 5 . Wave localization and soliton formation in this 
gap can be described within the framework of a discrete 
model of the tight-binding approximation, where the to- 
tal field is decomposed into a superposition of weakly 
overlapping modes <p n (x) of the individual waveguides 
in the form, £(x,z) = E n (z)cp n (x) ex.-p(i\ n z). Af- 
ter substituting this expression into Eq. QJ, we can ob- 
tain the modified discrete nonlinear Schrodinger (DNLS) 
equation for the normalized mode amplitudes E n , 



dz 



■E n + 1 )+ ln \E n \ 2 E n = 0. (2) 



Here A n characterizes the linear propagation constant of 
the mode guided by the n— th waveguide, j n are the ef- 
fective nonlinear coefficients, and A2 n +i = —p, ^2n = 0? 
where we assume appropriate normalization. 

Soliton solutions of Eq. (J2J) have the form E n = 
u n exp(i/3z), where (3 is the propagation constant, and 
u n is the soliton profile. Discrete bright solitons 
should have exponentially decaying tails, u±\ n + 2 m\ — 
u ±\n\ exp(— \2m\n) as n, m — > +oo. This condition is 
satisfied inside the band-gaps where Reft =^ 0. It can 
be demonstrated 5 - that k = (1/2 ) cosh -1 (— rj/2), where 
r] = 2 — /3(/3 + p). Then, discrete solitons can ap- 
pear in the semi-infinite total internal reflection gap, 
(3 > — (p/2) + yj (p/2) 2 + 4, and gap solitons can form 
inside the Bragg reflection gap for — p < (3 < 0. 

We can find an analytic estimate for the width of the 
gap soliton, a key parameter which can be controlled di- 
rectly in experiment. In the total internal reflection gap, 
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FIG. 1: Top: soliton power vs. propagation constant for A- 
type (dashed gray) and B-type (black) discrete gap solitons: 
solid — stable, dashed — unstable, and dotted — oscillatory 
unstable. Bottom: examples of the soliton profiles of (a) A 
and (b) B symmetries at (3 = —1; filled circles and squares 
mark the field amplitudes at the narrow and wide waveguides, 
respectively; p — 1.5. 



the width of the conventional discrete solitons can vanish 
for (3 — » +oo, until almost all the energy is confined in 
a single waveguide^. On the other hand, the situation 
is qualitatively different for discrete gap solitons, where 
the soliton width is bounded from below, and it can be 
estimated according to the linear tail asymptotic, 



A > 4/ max Re [«(/?)] 

/3,-p</3<0 



4/cosh _1 (l + p 2 / 8 )- (3) 



Thus, the minimum soliton width A is determined by the 
gap width (p), and the strongest localization is achieved 
at the middle of the gap ((3 ~ —p/2). For narrow gaps 
we have A > S/p ^> 1, and in this case the soliton dy- 
namics can be described by the coupled- mode equations 2 . 
However, for larger bandgaps the soliton can be localized 
at several waveguides (A ~ 1), and discreteness effects 
become very important. 

The soliton profiles and associated powers P = 
S n \ u n\ 2 are calculated numerically. We find that there 
exist two families of gap solitons, A and B, which bifur- 
cate from the lower gap edge ((3 = —p) in a self-focusing 
medium, see the example in Fig. Q As the power in- 
creases, the soliton width decreases, and it reaches the 
minimum value inside the gap, in agreement with the 
analytic estimate. Profiles of strongly localized discrete 
gap solitons are shown in Fig. ^bottom). 

A practically important question is soliton stability. 
We present the field in the form 
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and perform the linear stability analysis for the evolution 
of small-amplitude perturbations (v n and w n ) on top of 



FIG. 2: Example of (a) external and (b) internal resonances 
between an eigenmode of the linear eigenvalue problem and 
a bang-gap edges of the continuous spectrum that lead to an 
oscillatory instability of the odd-symmetry gap solitons. 



the soliton profile u n . After substituting Eq. (|2J) into 
Eq. J2J), we obtain a set of linear equations for v n and w n 
which possess localized solutions at discrete values of T. 
Solutions with ImT > indicate the soliton instability, 
as the perturbations grow exponentially. 

We find that the gap solitons of the A-type ex- 
hibit symmetry-breaking instabilities, and the instability 
growth rate becomes larger for strongly localized modes. 
This property is common for the systems with a broken 
translational symmetry, and it was demonstrated exper- 
imentally that such instabilities can strongly affect the 
soliton generation, and they can be utilized for the soli- 
ton steering 7 . 

Conventional discrete solitons of "odd" symmetry are 
always stable 7 . Similarly, we find that the B-type discrete 
gap solitons do not exhibit symmetry-breaking instabili- 
ties. However, our analysis reveals that these gap solitons 
become unstable above a critical power due to two types 
of the oscillatory instabilities. First, instability can ap- 
pear through an external resonance with the mode in a 
different gap, as illustrated in Fig. Efa). A similar in- 
stability mechanism was discussed earlier for nonlinear 
defect modes in a layered medium 4 . Second, the gap 
soliton can lose stability due to an internal resonance 
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FIG. 3: Excitation of discrete gap solitons by two interfering 
Gaussian beams: (a) self-focusing and soliton formation, (b) 
self-defocusing when the interference maxima are located at 
narrow and wide waveguides, respectively; (c) soliton instabil- 
ity at high powers; and (d) soliton steering due to power im- 
balance. Parameters are d\ — cfo — 2.5/xra, d s — 5/xra, 
refractive index no = 3.4, effective modulation Ano = 0.0014, 
vacuum wavelength 1.5/xra. Insert shows the generation ge- 
ometry. 



within the gap, see Fig. Efb). The latter scenario was 
first discovered for temporal Bragg solitons 8 . In both 
cases, instabilities are related to resonant excitation of 
the band modes, and the periodic beating is indicated by 
a non-zero real part of the instability eigenvalue T. 

Soliton excitation inside the frequency gaps of the fiber 
Bragg gratings is always accompanied by reflection of a 
substantial energy fraction of an incident pulse 2 . The 
same limitation is known for the side- on excitation of 
spatial gap solitons: an input beam is almost completely 
reflected when inclination angle is tuned deep inside the 
gap, and the soliton formation is observed only in the 
vicinity of gaps 3 . More efficient energy coupling into spa- 
tial gap solitons can be achieved by a head-on excitation, 
since the backward scattering is negligible. 

Spatial gap solitons can be generated by two Gaussian 
beams, which are tuned to the Bragg resonance and have 
opposite inclination angles 9 , as schematically illustrated 
in Fig. [3f insert). If the incident beams are identical, the 
input field can be written as 

Mx) = Ce-^-^l^ cos [tt(x - x s )/d] e ^+^V 2 . 



where x e is the beam center, d is width, C is amplitude, 
and (fi^ are the phases of interfering beams. Parameter 
%s = (<^2 — (pi)d/2ir defines a shift of the interference 
pattern due to the relative phase between the beams. 

Our analysis of the gap soliton solutions revealed that 
the soliton energy is mainly localized at the narrow 
waveguides, see Fig. ^bottom). We perform numeri- 
cal simulations using Eq. and confirm that, when 
the interference maxima are at the narrow waveguides, 
the two Gaussian beams experience mutual focusing and 
a gap soliton forms, see Fig. Eta). The optimal input 
powers and beam widths can be chosen to reduce sig- 
nificantly the emission of radiation. On the other hand, 
if the interference maxima appear at wide waveguides, 
nonlinear beam interaction accelerates the beam spread- 
ing and break-up [Fig. EDO]- The soliton generation 
by two beams becomes impossible for higher powers [cf. 
Figs. Et a ) andEt c )]? an d the generated beam exhibits 
periodic beating during its propagation along the array. 
This happens due to the development of soliton insta- 
bility which limits mutual beam focusing through a reso- 
nant excitation of modes in different bands, in agreement 
with our linear stability analysis. Finally, we illustrate 
how the soliton motion can be induced by varying the 
power imbalance of the input beams [Fig.Etd)]. 

In conclusion, we have analyzed the basic properties of 
discrete gap solitons using engineered binary arrays as an 
example. For the first time to our knowledge, we have de- 
scribed the families of discrete gap solitons and analyzed 
their stability. We have revealed two basic mechanisms of 
the soliton instability and discussed its connection with 
mutual beam focusing and gap soliton generation in pe- 
riodic photonic structures. 

The work was supported by the Australian Research 
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